Abstract 



In the selective withdrawal experiment fluid is withdrawn through a 
tube with its tip suspended a distance S above an unperturbed two-fluid 
interface. At low withdrawal rates, Q, the interface forms a steady state 
hump and only the upper fluid is withdrawn. When Q is increased (or 
S decreased), the interface undergoes a transition so that the lower fluid 
is entrained with the upper one, forming a thin steady-state spout. Near 
this discontinuous transition the hump curvature becomes very large and 
displays power-law scaling behavior. This scaling is used to show that 
steady-state profiles for humps at different flow rates and tube heights 
can all be scaled onto a single similarity profile. 
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Is it possible to classify topological transitions in nonlinear fluid systems 
I II in the same manner as one classifies thermodynamic transitions? When 
the topological transition involves formation of a singularity in the fluid flows 
or interface shapes, a similarity solution can provide a simplified description 
of the flows and help make such a classification Q. A crucial component to 
this approach involves determining how characteristic physical quantities and 
lengths describing the fluid system scale near the singularity. In many cases 
these singularities manifest themselves in the transition dynamics Q, H, ^ 
and do not appear in the steady state flows. Here, we report on steady-state 
interface profiles near the topological transition associated with the selective 
withdrawal experiment. Despite the transition being discontinuous, scaling of 
the interface is observed as the transition is approached. 

In the selective withdrawal experiment a tube is immersed in a filled con- 
tainer so that its tip is suspended a height S above an unperturbed interface 
separating two immiscible fluids. When fluid is pumped out through the tube 
at low flow rates, Q, only the upper fluid is withdrawn and the interface is 
deformed into an axi-symmetric steady-state hump (Fig. |]) due to the flows 
in the upper fluid. The hump grows in height and curvature as Q increases or 
S decreases until the flows undergo a transition where the lower fluid becomes 
entrained in a thin axi-symmetric spout along with the upper fluid. The two- 
fluid interface becomes unbounded in the vertical direction thus changing the 
topology of the steady state. Once the spout has formed, an increase in Q or 
decrease in S thickens the spout. 

The interfacial profiles at different flow rates and tube heights are recorded. 
Near the transition, the steady-state radius of curvature of the hump tip is 
orders of magnitude smaller than the length scales characterizing the boundary 
conditions (e.g. the tube diameter, D). This separation of length scales suggests 
that a similarity analysis of the steady-state hump profiles might be possible. 
However, for the range of parameters explored thus far, even when the system 
is arbitrarily close to the transition from hump to spout, the mean curvature 
of the hump tip, k, while large, remains finite. Nevertheless, by fixing S and 
looking at the steady-state profiles as Q is increased, we observe that both the 
hump curvature and height display scaling behavior characteristic of systems 
approaching a singularity. Since the divergence is cut off before a singularity is 
reached this transition appears to be "weakly-first-order." 
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Figure 1: Photograph of two- fluid interface forming a hump. The parameters for 
this problem include the height of the tube from the unperturbed interface, S, 
withdrawal rate, Q, surface tension of the two-fluid interface, 7, fluid densities, 
Pu and pi, fluid viscosities, rju and r]i, tube diameter, D, fluid height above 
the interface, fluid height below the interface, container size, and surfactant 
concentration. Here we restrict the investigation to the S vs. Q parameter 
space and defer discussion of the remaining parameters to [nOj . 
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As shown in Fig. the parameters important for this experiment are the 
upper and lower fluid viscosities and densities (ryu, 77;, p„, p/), interfacial tension 
(7), orifice diameter (£>), tube height {S), and flow rate (Q). In looking for 
scaling of the steady-state profiles, care must be taken to design an experimental 
apparatus capable of isolating the profiles near the transition. Experiments were 
performed in large tanks (30cm x 30cm x 30cm) capable of holding fluid layers 
that were each about 12 cm in height. To ensure that the upper fluid level 
remained constant, the withdrawn fluid was recycled back into the container 
(the bottom fluid layer thickness remains constant when the system is in the 
hump state). Steady withdrawal rates were achieved by using a gear pump. 
We verified that for the tube diameter {D = 0.16 cm), tube heights (0.1 cm 
< < 1.1 cm), and flow rates (Q < 10 ml/sec) used in the experiments, the 
container walls were sufficiently distant and the fiuid layers sufficiently thick so 
as not to affect the fiows [0. We measured the upper(Heavy Mineral Oil) and 
lower (Glycerin- Water mixture) fluid viscosities to be rju = 2.29 St and rji — 1.90 
St, the upper and lower fluid densities to be p„ ~ 0.88 g/ml and pi = 1.24 g/ml, 
and the surface tension |l^, |l^ to be 7 = 31 dynes/cm. Attempts to increase 
K near the transition by decreasing 7 cause fluid mixing and result in a diffuse 
interface at high shear rates so that some fraction of the lower fluid is always 
being withdrawn ]l3[ |. 

While many of the parameters mentioned influence the flows, our under- 
standing of the scaling behavior can be conveyed by focusing on S and Q. We 
can fix Q and track the development of the hump profiles as a function of S. 
Below the tube height, S'u, the hump is unstable and undergoes a transition to 
a spout. Figure | shows that 5„ oc qO-30±o.05 q ^j^g transition is 

hysteretic: the value of S where the spout becomes unstable and decays back 
into the hump is larger than Su- We define the difference of the two heights 
or hysteresis as AS". Figure |^ indicates that the data is consistent with an ex- 
ponential decrease: AS" = 0.04 exp"'^/^ ^^^. For Q > 0.1 ml/sec, AS" was too 
small to measure. 

In order to measure the mean curvature of the hump tip, k, we first fit the 
tip of the recorded profile with a Gaussian and then calculate the curvature 
of the fitting function at the hump tip. Figure || also shows the hump height, 
hu, and mean radius of curvature, l/nu, at the transition as a function of Q. 
The dramatic decrease in AS" coincides with the onset of a flat asymptotic 
dependence for I/ku at Q > 0.1 ml/sec. We quantify this correlation, by fitting 
the curvature data with the form 1/ku = 0.02 -I- 0.32 exp~*5/° °^^ which has the 
same exponential decay with Q as does AS'. For Q > 0.1, ml/sec we find both 
Ku, and hu to be independent of the orifice diameter D ||l^. We restrict our 
scaling analysis to this regime. Jisf 

Figure 3a plots k (where < k < k^), versus Q for six data sets correspond- 
ing to different values of S. As shown in the inset of Fig. 3a, all fifteen data sets 
display a power-law divergence for k, as Q approaches Qc (a fitting parameter). 
While the power-law exponents remain constant as S is varied, the prefactors 
to the power laws, Ck{S), vary slightly with S and are scaled out in the inset. 
Figure 3b plots the hump height, hmax, versus Q. The inset to Fig. 3b shows 
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Figure 2: Plots of the transition tube height, Su, hysteresis, AS', hump height, 
hu, and mean radius of curvature 1/ku, as a function of Q. We find that 
Su oc qo.30±o.05 ^gQijfj line). The AS" data is fit with an exponential decay 
(dashed) of the form: AS = 0.04exp~'5/° °^^ although the range is not sufficient 
to exclude a power-law decay. We fit (dash dot), with the form: = 

0.02 + 0.32 exp-'3/°°32. 
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that as Q approaches Qc (obtained from Fig. 3a inset), the hump height ap- 
proaches he (a fitting parameter) as a power law. Once again, the power-law 
exponents for these data sets are the same over this range of S. The prefactors, 
Ch{S), are scaled out in the inset. Note that Qc changes with S indicating that 
the system can approach a continuous line of divergences. Combining the two 
scaling dependencies in Fig. 3c, we plot {he — hmax) / hmax versus the normalized 
curvature, k/ti. We find that {he — hmax)/hmax scales as (K/n)°-^^^°"^ indi- 
cating that even though both he and the power-law prefactor, n, change with 
S, the power-law exponents are independent of S for this range of tube heights. 
Note that n{S) = Ch{S)[ci^{S)°-^^]. The transition cuts off the evolution of the 
hump states making it impossible for the system to approach arbitrarily close 
to the singularity and limiting the precision with which we can determine the 
exponents. 

The scaling observed for hmax and k suggests that the hump profiles should 
display universal behavior as hmax nears he- The quantities l/(/«/n) and {he — 
hmax) /hmax track how quickly the radial and axial length scales decrease as the 
system approaches the singularity and are therefore used to scale the profiles. 
We define the scaled variables: 

H{R) ^ ^ ' and i? = — , (1) 

f^c f^max ^ 

where h{r) is the hump profile and he is taken from Fig. 3. The transformation 
shifts the profiles so that under scaling the singularity occurs at the origin and 
the maximum hump heights occurs at = 1 and R = 0. Figure 4 shows eight 
scaled profiles for the S — 0.830 cm data set. In the bottom inset we overlay the 
hump profiles scaled in the main figure. We find excellent collapse of the profiles. 
The solid line in Fig. 4 is a power law that fits the data in the region beyond 
the parabolic tip. The picture that emerges is of a parabolic tip region which 
decreases in its radial scale and is simultaneously pulled towards the singularity 
in the axial direction leaving behind it a power-law profile with an exponent 
of 0.72 ± 0.08. This exponent is within error (although slightly smaller) of the 
exponent observed in the scaling relation of Fig. 3c which predicts an exponent 
of 0.85 ±0.09. 

Typically, the observed scaling dependencies in these types of problems result 
from the local stress balance. A scaling analysis where the viscous stresses 
of the upper and lower fluids balance the stress arising from the interfacial 
curvature predicts linear scaling dependencies and conical profile shapes. The 
non-linearity of the observed dependencies indicates that either a different stress 
balance governs the flows (e.g. only viscous stress due to upper fluid balances 
stress due to the interface curvature) or that non-local effects are coupling into 
the solution. A more detailed discussion can be found in | pO[ . 

Finally, we compare the similarity curves for five different tube heights in 
the upper right inset of Fig. 4. The profiles all display the same power-law 
dependence. Within error, the normalized curvature k/h (taken from Fig. 3c) 
can be used to scale the radial components of these profiles and obtain good 
collapse. In Fig. 3c we find that the normalization prefactors, n, decrease as 
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Figure 3: Scaling for the mean hump tip curvature, k, and height hmax- 3a 
plots K, vs. Q for six tube heights. Each curve displays diverging behavior with 
increasing Q. For each tube height we choose a critical flow rate Qc as a fitting 
parameter and show (3a inset) that as Q approaches Qc the curvatures increase 
as power laws. The prefactors to the curvature power laws, Ck{S), are scaled 
out in the inset. Figure 3b plots h^ax vs. Q for six tube heights. For each S 
we choose a critical hump height he as a fitting parameter. As Q approaches Qc 
(taken from 3a inset), the hump heights approach the critical heights as power- 
laws. The prefactors, Ch{S), to the power laws in the inset are scaled out. 
Figure 3c plots {he — hmax) /hmax vs. k/u for the entire data set corresponding 
to fifteen different straw heights. The prefactor n roughly decreases as exp"^-^'^. 
The line corresponds to a power-law with an exponent of -0.85. 
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exp~^-^'^. Here, we correlate this decrease with the observation that the profiles 
become shallower at larger S. The points of deviation for the S = 0.255 cm and 
0.381 cm profiles mark the transition from the similarity regime to the matching 
regime beyond which the profiles become asymptotically flat. At large enough 
radii all of the scaled profiles display these deviations. 

We have shown that in the Q > 0.1, ml/sec regime, a similarity analysis 
can be used to describe the flows near the selective- withdrawal transition p6{ . 
We have observed power-law scaling of the hump height and curvature (Fig. 
3) and used these scaling relations to collapse the hump profiles at different 
flow rates and tube heights onto a single universal curve (Fig. 4). However, 
the origin of the saturation of k at large Q remains an important unexplained 
problem. Further insight into this cutoff behavior may be gained by comparing 
with an analogous two-dimensional (2-D) problem which roughly corresponds 
to replacing the tube with a line sink. Jeong and MofFatt showed that in 
an idealized case where the bottom fluid is inviscid while the top fluid is very 
viscous, the 2-D hump interface forms a 2-D cusp singularity above sufficiently 
high withdrawal rates. Recently, Eggers showed that the solution changes 
when the lower fluid has a finite viscosity; the system no longer manifests a 
singularity jl^. Instead, the approach to the singularity is cut off and the 
system undergoes a transition to a different steady state. In this new state, a 
sheet of the lower fluid is entrained along with the upper fluid into the line sink. 
However, the finite lower fluid viscosity prevents the hump profiles from scaling 
onto a similarity solution. 

Here, we have shown that for the three-dimensional selective withdrawal sys- 
tem even when both fluids are viscous (77 > 1 St for both fluids) the effects of 
a singularity manifest themselves in the scaling of the hump profiles. Further- 
more, preliminary experiments show that a reduction of the lower fluid viscosity 
to 0.01 St has little effect in determining the flnal curvature of the hump tip 
or equivalently, how close the system is to forming a cusp. This suggests that 
for our 3-D problem, either the effects of the lower fluid viscosity enter as a 
higher-order perturbation to the proflle shapes, or a different mechanism un- 
derlies the avoidance of the cusp formation. If the latter scenario is correct, it 
may be possible for the system to manifest the singularity at a flnite lower fluid 
viscosity. In either case, determining which variables affect how close the sys- 
tem approaches the singularity would allow for control of the maximum hump 
curvature and minimum spout diameter. This control could then be used to ad- 
vance technologies such as coating microparticles pof , creating mono-dispersed 



micro-spheres [£1[, and emulsification through tip streaming [y_8, 22 which take 
advantage of the selective withdrawal geometry. 
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Figure 4: The scaled hump profiles. The lower inset shows eight profiles taken 
from the S = 0.830 data set. The main figure shows the same profiles after 
scaling. The solid line corresponds to a power law of the form i?"-^^. In the 
upper inset we compare the universal curves for the S = 0.830 cm, 0.613 cm, 
0.508 cm, 0.381 cm, 0.255 cm data sets. 
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